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Abstract: 

We show the equivalence of a previously derived exact expression for the response of a 
non-relativistic system with harmonic forces and an infinite sum of weighted J-functions cor- 
responding to the spectrum. We forward arguments, indicating that the Gersch-Rodriguez- 
Smith 1/5-expansion of the response does not converge and prove that this expansion is an 
asymptotic series. 
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In recent years there has been renewed interest in work by Gersch, Rodriguez and Smith 
(GRS). Those authors formally derived for large momentum transfer q, a. 1/q expansion of 
the dominant incoherent part of the total non-relativistic (NR) response, S{q,u)), which is 
non-perturbative in V^. This approach has recently been applied, even if the underlying 
forces are singular and attractive (i.e. confining) ^. The incoherent part of the total reduced 
response (l){q,y) = {q/m)S{q,(jj) = + (j)''°^{q,y) may thus be expressed as 

oo 

^i^^o'^{q,y) = Y^{m/q)''Fk{y), (1) 

where the energy loss w has been replaced by an alternative variable 

y = -q/2 + mw/q (2) 

Up to now a detailed analysis of the convergence of GRS expansion of the response for NR 
systems, interacting either by regular or confining forces, is lacking. In the present note 
we consider the possibihty that, even if this series does not converge, a finite number of 
terms may still describe the response at high q to any desired accuracy, as is the case for an 
asymptotic series. 

Without reference to specific NR dynamics we consider the following three equivalent 
forms for the response of a target of mass tua, composed of A constituents with equal mass 
m 

Siq,ij) = A-^ Yl \^On{q)\''SiiJ - q^ l^rriA - Eno) (3a) 

n 

= {TrA)-\0\pllmiuj + Eo-T-V-iv)-' p,|0) (36) 

/OO 
dse'^y e'"«/2(0|pt exp 
-oo 

where J^on{q) = (0|pj|n) is the inelastic form factor between the ground state |0) and excited 
states |n), Eno are intrinsic excitation energies without the target recoil energy q'^/2MA, which 
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{Eo - H) 



(3c) 



appears separately in the energy conserving 5 function in (3a), and pq is the Fourier transform 
of the charge density. A formal summation over the excited states in (3a) produces (3b) in 
terms of the Green's function with T, V the total kinetic and potential energy. Eq. (3c) is 
the Fourier transform of (3b) and H = T + V the total Hamiltonian of the system. Clearly it 
is meaningful only in the distribution theory. 

We now concentrate on systems with confining interactions and shall be led to questions 
of equivalence of the different expressions (3). We shall illustrate our reasoning on the 
example of harmonic confining forces, for which an exact solution for the response in terms 
of a single integration has recently been derived ^. It suffices to treat the case of two particles 
in a HO well. With /3 = (mf2/2)^/2 ^j^g inverse oscillator length and a = 2f3'^s/q one finds for 
the reduced response 

^i--oh^g^y^ ^ j _exp(zys)$(a,s) (4a) 

^co.(,, = £ I exp (^) ] -P {^^^^4 + l) } (4^) 

with 

$(a,«)=exp|-^^-^j f^P|-^^(,— -Vl 

The exponents in Eq. (5) are readily seen to be sums of power series in a and thus permit 
an expansion of ^{a, s) in series of 1/q, except for g' = 

oo 

^a,s) = 2 7Tj2a/Q)''hk{s) (6) 

When permitted, interchange of summation and term by term integration leads to the GRS 
sum (1). The HO example is particularly useful for the discussion of two points: 

1) With a priori knowledge of the equivalence of the expressions (3), can one directly 
demonstrate that the integrands in Eqs. (4), analytic except for q = 0, lead to a weighted 
sum of ^-functions? 
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2) The GRS expansion is formally generated through term by term integration of Eq. 
(4a) and it consists of a series of regular functions, (cf. Ref. 3, Eq. (3.1)). The sequence of 
the partial sums of this expansion does not seem to tend to a series of equally spaced peaks, 
required to obtain in the limit Eq. (3a), i.e. a sum of weighted, equally spaced 5- functions. 
In view of this unlikely equivalence, what is the nature of the GRS series (1)? 

We start with the first point and rewrite (4a) 

r°° fin 

r^^^q, y) = exp(-gV8/3') j ^ exp(zys) exp(zgs/4) D{q, s) 
D{q, s) = exp 

The function D{q,s) is periodic in 2p'^s/q and thus permits the Fourier expansion 

+ 0O 

D{q,s)= J2 dniq,P)e-^'^f''^/^ (8) 

n= — oo 

with expansion coefficients 
and 

= (-1)" J dz ^-("+i)e-^ = ^ . (10) 

The integration in Eq. (10) is performed on a circle centered in the origin and of radius 
|c| = g^/S/J^. Eq. (9) can be readily obtained from Eq. (10) with the substitutions z = 
ce-*" (27r > a > 0) and c = -q'^/8/3^. 

In order to show the equivalence of Eqs. (3a) and (3c), we recall that Eqs. (3c), (4a) 
and (7) have meaning only in the framework of distributions. There we need the notion of 
'good' functions g{y), i.e. functions which are everywhere differentiable any number of times, 
and which together with their derivatives vanish as \y\ oo faster than any power of l/\y\ 
For any such g, Eq. (7) implies 

cl)'"''°\q, y) 9{y) dy = (27r)-i exp(-gV8/32) / e^^/^ D{q, s) g{s) ds, (11) 

-OO J —OO 



where g{s) is the Fourier transform of g{y). We note moreover 

N 
n=-N 

and that the right hand side of Eq. (12) is independent of N. The Lebesgue dominated 
convergence theorem '^^ then permits inversion of the order of integration and summation of 
the series in Eq. (7) and one readily finds 

5{ui - q^/Am - nCl) 

(13) 

5[y + q/4-{2p^n/q)] 

A similar expression can be derived for the coherent part and one checks that the weights 
of the 5- functions are the exact squared inelastic form factors for the HO^. Notice that for 
w > 0, only n > 0, i.e. the actual HO spectrum, contributes to (13). 

In spite of its utter simplicity, the above explicit demonstration of equivalence is not 
trivial at all, nor does it serve as a model for other confining interactions. Only for a very 
limited number of Hamiltonians can one give closed forms for the exact response, using 
either exact Green's function^ or operator techniques^. In addition, for no other interaction 
can the corresponding function D{q, s) in (8) be periodic: it would lead to an equally spaced 
spectrum, characteristic for the HO. In general a regular Hamiltonian has a finite number 
of unequally spaced discrete eigenvalues. In any case the proof above for the equivalence of 
Eqs.(3a) and (3c) is not readily seen to be of a general type. 

We return to Eq. (4a). We already observed that the integrand in Eq. (4a) permits 
a 1/g-expansion, which produces the GRS series if the order of summation and integration 
may be interchanged. However, we have not been able to prove a condition like Eq. (12) 
for the partial sum of the power expansion of $(a, s), Eq. (6). A finite number of terms in 
the 1/q expansion may always be integrated term by term, thus leading to a truncated GRS 
series and this has been standard practice for systems with regular interactions. A prime 
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8/3^ 



(12) 



^ ' q ^ n\ 



q 

8/32 



example is the successful treatment of liquid ^He . However, the crucial question is whether 
the above interchange is permitted for the full series. If the GRS series does not converge, 
it may still be an asymptotic one ^ and we now show that this is the case for the HO. 
Consider the sequence {Sn{q,y)} 



n p+oc 

Sn{q,y) = / hk{s)e'yUs. 

fe=o 



(14) 



All integrals in Eq. (14) are well defined, since hk{s) is a linear combination of powers of s 
times exp(-/?^s2/4) and therefore hk{s) e jC^. We wish to show that the sequence {Sn{q,y)} 
defines an asymptotic series for To this end we have to show that for any fixed 

n and any good function g{y), for any given e > a Q > exists such that, ior q> Q one has 

An{q) < e '^'^ with 

/ + 00 /■ + 00 

/ Sn{q,y)g{y)dy 
-oo J — oo 

/ + 00 /' + 00 

Y,iyq)''hk{s)g{s)ds - / Y.{l/q)''hk{s)g{s)ds 
-°° k=0 k=0 

In the second term of the last step of Eq. (15) the order of integrations has been exchanged 
by using the Fubini theorem^*". Eq. (15) can be put in the following form 



An{q) =q" 



= q" 



(15) 



K{q) 



'i>^''+'\ao,s)g{s)ds 



(16) 



. (n + 1)! 

where $("+i)(a,s) = (5/aa:)("+i)$(a, s) , x = l/q = a/{20'^s), and = '^P'^s/qo with Qo > q. If one 
can show that for any a e (-oo, +oo) 



$("+i)(a,s) 



(n+l) 



is), 



(17) 



with 'P(„+i)(s) diverging at most as a power of s, our theorem is demonstrated. It is easily 
seen that the derivative $("+^)(a, .s) is a product of $(q:,,s) and combinations of £-th. order 
derivatives ( £ < n+l) of a{a) = sin(a/2)/(a/2) and b{a) = (sina/a- l)/a with respect to a, with 
coefficients depending on powers of s with bounded exponents for a given n. 
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We thus have to study the behaviour of a(a), b{a) and their derivatives for a g (-00, +00). 
Consider separately the two sets \a\ < 1 and \a\ > 1. In the range [-1, 1] one may expand 
a{a) and 6(a) as a power series in a, which together with the series defining their derivatives, 
converge uniformly in the closed set [-1, 1] and their absolute values have there finite maxima. 
For \a\ > 1, the functions a{a), b{a) and their derivatives can be expressed as combinations of 
sina, cosa and powers of l/\a\ < 1. One can therefore find quantities, independent of a, which 
are larger than the absolute values of a{a), 6(a) and their derivatives for any a e {-00, +00). 
Consequently (17) holds and our theorem is demonstrated. 

We summarize: Although the GRS expansion for the response of NR systems turns out 
to be a useful tool, we are not aware of an actual proof that the above expansion for systems, 
interacting through regular or confining forces is a convergent, or at least an asymptotic one. 
In this note we have proved that the latter is the case for the 1/g-expansion of the response 
for an harmonic oscillator. 
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